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1. INTRODUCTION 

A model of high frequency radio (HF) skywave modes under ambient 
conditions is basic to the study of HF system response to disturbed condi¬ 
tions, such as those which would be encountered in a nuclear conflict. 

Traditionally, HF skywave vulnerability has been calculated under 
the assumptions that only D-layer absorption is important, that ambient 
propagating modes are an adequate approximation to modes which propagate 
after nuclear explosions, and that a careful description of late time (ten 
minutes to hours) debris motion is unnecessary so long as an assured defense 
conservative model is used. 

Recent interest in adaptive HF systems for the Air Force and 
MEECN, as well as interest in the closely related class of meteor burst 
communication systems, has forced a re-evaluation of all three of the above 
assumptions. More realistic models of late time nuclear debris motion (the 
source of nuclear radiations which ultimately cause D-layer absorption) 
have recently been developed 1 ’ 2 to allow improved absorption calculations. 

In this report we describe a method of ambient HF skywave mode calculation 
designed to readily accept disturbed ionospheric conditions, and thus to 
allow inclusion of some nuclear disturbances important to propagation 
geometry, necessary for more accurate calculation of D-layer absorption 
effects as well as for removal of the second and third traditional assump¬ 
tions quoted above. 


3 






1 


The necessity for development of an ambient skywave propagation 
mode calculation requires justification. Elaborate techniques for detailed 
calculation exist in the scientific community, 3 which have been adapted 
into a nuclear effects code. 1 * Less elaborate techniques exist for long 
term average mode prediction, 5 to which simplified fits have been made, 
and these fits are commonly used in the nuclear effects community. 6 To 
understand the requirement for a technique of intermediate complexity re¬ 
quires a brief discussion of anticipated nuclear effects of HF propaga¬ 
tion in relation to characteristics of available calculational techniques. 

Anticipated nuclear effects which may conveniently be taken 
into account by alteration of normal ionospheric parameters are changes 
in ambient layer (i.e., F2, FI, or E-layer) electron density, thickness, 
height, and tilt due to: (1) propagating acoustic gravity waves, (2) ducted 
waves, (3) burst region disturbance (F-layer depletion, etc.), or (4) bomb 
or debris irradiation. Another set of effects exists (M-modes, fireball 
modes, plume modes, etc.) which are more conveniently calculated indepen¬ 
dently and which will not concern us here. To calculate, even crudely, 
effects of the above list on propagation geometry, thus on existence and 
efficacy of ambient-like modes and creation of new modes, the ambient 
mode calculation must readily take into account ionospheric variations all 
along the propagation path in terms of variability of reflection height 
and ionospheric tilt, both along and across the great circle path. In 
fact, non great circle propagation must be treated. The current ITS78 
code, 5 the community standard, evaluates ionospheric conditions only at the 
great circle path center, derives from that a single reflection height, 
and does not consider tilt or off-great circle propagation. This is satis¬ 
factory for long term average predictions, the objective of that code, but 
does not lend itself to situations similar to nuclear effects where a portion 
of the path is strongly disturbed. The WRECS code, 6 using a fit to ITS78 
predictions, allows for two values of reflection height corresponding to 
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day or night but is otherwise restricted to ITS78 predictions. The NUCOM 
code 4 uses ray tracing techniques and thus can be much more versatile, but 
such techniques are expensive for our purposes. 


The method developed here to meet requirements of nuclear effects 
studies is approximate but more accurate than justified by predictability 
of either the ambient or nuclear disturbed ionosphere. It is more expen¬ 
sive in computer time than one would prefer but still vastly cheaper than 
normal ray trace techniques; sophisticated in that it accounts for ionos¬ 
pheric variation, tilt, and off-great circle propagation; but in its pres¬ 
ent implementation rudimentary in that only two modes can be calculated for 
a given number of hops, no extraordinary ray modes are calculated, etc. The 
restriction to two modes per hop number could be readily remedied and this, 
at least, is hoped to be done soon. The basic technique is believed to be 
sufficiently versatile that extension toward more realistic treatment is 
feasible. As reported here the method has been develoepd to a useful state, 
one that is more realistic than others in many important respects, as well as 
one that is reasonably inexpensive in application. 

An outline of the overall calculation is given in the next sec¬ 
tion and specific components discussed in following sections. 






OUTLINE OF THE CALCULATION 

The calculation is naturally broken into two logical modules, 
each of which consists of several component subroutines. 

2.1 Ionospheric Data 

The first module is exercised once only, before any propagation 
calculations are attempted. It uses an ionospheric model to create a data 
table of ionospheric parameters for each desired transmitter-receiver link. 
Each such data table consists of three sets (one for each ionospheric layer, 
E, FI, F2) of three parameters describing parabolic ionospheric layers 
above several points on the great circle path at a number of times. The 
three layer parameters are: (1) maximum electron density, (2) thickness 
of the layer, and (3) altitude of the maximum electron density. The 
spatial points are arbitrary in number and spaced equally along the great 
circle path. The number of times is also arbitrary and determines a number 
of equal segments of a twenty-four hour period at which each layer will be 
specified at each point along the great circle. Normally we choose twenty- 
four times, one for each hour of the day (plus one repeated entry to allow 
the table to begin at 0 hours and end at 24 hours), and seven space points; 
which yield six equally spaced points along the path (both end points are 
included), normally these great circle points are separated by less than 
15° of arc each (that is, corresponding to less than an hour apart local 
time). With ionospheric models at our disposal 5,7 one hour equivalent seg¬ 
ments have been found to be about twice as dense as necessary to allow 
interpolation in both space and time to the limit of model fidelity. 

The purpose of this array of (3 layers) * (3 parameters) x 
(typically 7 space points) x (typically 25 time points) = 1575 numbers 
per link is to allow rapid interpolation in space and time for ionospheric 
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quantities without having to repeatedly calculate potent! ->lly complex ionos¬ 
pheric models for identical or similar conditions. Heights and thicknesses 
are interpolated linearly between adjacent points and electron density 
logarithmically. In Section 5 we shall show that despite appearances, 
this form of ionospheric data table does not restrict consideration to 
great circle paths. 

A typical comparison of the three parabolic layer ionospheric 
approximation used here to the ITS78 model, 5 the Aerospace Model 7 and to data 7 
is shown in Figure 1. Cases can be found in which the three-layer approxima¬ 
tion fits the model somewhat better and somewhat worse than Figure 1. It is 
also easy to find both better and worse comparisons between the models and 
data. Figure 1 illustrates the general result, which appears to be valid 
for all ionosphere models which might be used as a basis for the fit, that 
error resulting from the three parabolic layer approximation is small com¬ 
pared to error due to the uncertain behavior of the real ionosphere. 

2.2 Mode Calculation 

The second module calculates mode geometry. Inputs to this module 
are a specified time, ray type, link, frequency, and number of hops. All of 
these quantities are arbitrary except for obvious restrictions, frequency 
and number of hops must be greater than zero; and except for a bow to the 
nature of Fortran compilers, presently we dimension for no more than eight 
hops. The ray type specifies either the low ray or high (Pederson) ray. 
Outputs of this module are either specification of mode geometry or infor¬ 
mation concerning the reason no such mode exists. 

Mode geometry is specified in terms of the ends of a number of 
straight line segments (passes) connecting reflection points. A one-hop 
mode is specified by latitude, longitude, and altitude above the earth's 
surface of three line segment ends; (1) the transmitter coordinates. 
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Altitude, km 



Figure 1. Typical ambient ionosphere. The HFNET three parabolic layer fit 
is compared to the Aerospace model from which it is derived, to 
the ITS78 model, and to data for 1200 local time above Arecibo, 
Puerto Rico, Dec. 1968. Smoothed Zurich sunspot number is 110. 







(2) the equivalent triangular ionospheric reflection point, and (3) the 
receiver coordinates. A two-hop mode requires specification of coordi¬ 
nates for five points; the two end points, two ionospheric reflection 
points, and an intermediate surface reflection point. In general three- 
space coordinates must be specified for 2n + 1 points for an n-hop mode. 

A mode fails to exist if: (1) the F2 layer fails to reflect 
the ray on any hop, (2) transmitted launch angle for the mode is below 
a specified value, (3) received angle is below a specified value, 

(4) ionospheric tilt causes the ray to be reflected back upon itself, or 

(5) ionospheric tilt causes the ray to be deflected forward at too great 
an angle. 


Steps in calculation of n-hop mode geometry are as follows: 

1. The great circle path is divided into n equal segments, 
one for each hop. 

2. For each segment a one-hop mode is calculated according 
to local ionospheric conditions as follows: 

a. Assume reflection occurs at the F2-layer; calculate 
the equivalent triangular reflection altitude, h r , 
for a parabolic layer by iteration patterned after 
Reference 5. 

b. If successful, calculate refraction effect due to 
the FI layer. If the FI layer reflects the ray, 
carry out (a) for the FI layer. 

c. Calculate refraction effect due to the E-layer. If 
the E-layer reflects the ray, carry out (a) for 
the E-layer. 


9 







The situation at the end of step (2) is depicted in the sketch, 
provided all hops were successfully reflected. At this point in the cal¬ 
culation the mode is specified by latitude and longitude end points equally 
spaced along the great circle path but each ionospheric reflection height 
is different, a violation of the laws of optics. 

3. Adjust the pass ranges, r^, so that incident and reflected 
angles are equal, accounting for differences in h^. 

4. Calculate ionospheric tilt angles both along the great 
circle and across the great circle at each ionospheric 
reflection point. 

5. Adjust all pass coordinates (latitude, longitude, altitude) 
to account for tilt along the great circle to first order 
in tilt angle (that is, in the vicinity of each reflection 
point, the tilt is constant). If the ray is reflected back 
upon itself at any point, or if the cumulative effect of 
all tilts cause the ray to land too far beyond the re¬ 
ceiver, reject the mode. 
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6. Calculate off-great circle displacements of all pass coordi¬ 
nates and new reflection altitudes due to across great circle 
tilt. 

7. Find launch and received ray elevation relative to the horizon, 
and azimuth relative to the great circle path. If either 
elevation is below a specified limit, reject the mode. 

At the end of this process an approximation has been found to a 
"good" propagating mode with specified number of hops, within the limitations 
indicated above. 

Possible paths which include hops of the form indicated in the 
sketch are not found. The situation depicted in (a) is rejected because 




r 2 is negative, the ray was reflected back upon itself. Within the mathe¬ 
matical restrictions of first order tilt angle treatment this situation can 
only arise when h r is greater than the chord between the two ground points, 
and can only represent a reasonably plausible situation when the second point 
represents receiver or transmitter. These conditions are too severe for am¬ 
bient skywave propagation. Situation (b) is rejected because the calculation 
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does not allow for existence of the second reflection point, 2, and assumes 
that the ray continues along the dotted path. This situation seems unlikely 
under ambient conditions but could occur in a disturbed environment. It 
must be handled, if at all, in future as a "disturbed" mode, with appro¬ 
priate logic specified. 

The reader will note that no step between (1) and (7) was iterated. 
The only iterations occur in finding single-hop solutions to reflection from 
a single parabolic layer in step (2). With this exception the mathematical 
treatment runs "open loop" and could be made more rigorous by iterating on 
each refraction correction, after range correction for differing reflection 
heights, and after corrections for tilt. Such a procedure would be pro¬ 
hibitively expensive in calculation time and would produce numerical results 
significantly different from the present technique only for cases where some 
reflection occurred very near a critical point (e.g., almost or barely re¬ 
flected off some layer). Predictability of either natural or disturbed 
ionospheric conditions is not sufficiently accurate to justify such fine 
tuning. 


Of much greater importance is the fact that this method does not 
calculate all possible "sensible" rays. This defect is illustrated in the 
sketch, in plane geometry for illustrative purposes. 



Ground 
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Our technique first finds a tentative (low or Pederson) ray which 
reflects off the highest available layer, the dashed ray of the sketch which 
reflects in the vicinity of point A. Next the refraction angle due to pas¬ 
sage through the lower layer at point B is calculated. If refraction is 
not large enough to reflect the ray at B, a correction to the effective 
triangular reflection point A is introduced and calculation proceeds to the 
next step. At this point our technique has irretrievably missed a potential 
lower ray which reflects at C. Under almost all ambient conditions the lowest 
ray suffers greater absorption and also is a more likely candidate for rejec¬ 
tion due to low launch angle, thus is not likely to be useful or to cause 
multipath problems. It is unusual but possible that nuclear effects could 
preferentially interfere with the upper ray enough to make the ray at C the 
best mode. This is a weakness in current logic which probably should be fixed. 

On the other hand, if the ray is reflected at B then the tech¬ 
nique shifts attention to the lower layer and calculates only the ray 
which reflects at C. Another option would be to search at higher angles 
for a ray which is strongly refracted at D, closer to the critical angle 
for the frequency if such exists, and which traverses the path indicated 
by the light solid line, a high ray. 

Thus the technique will select preferentially the dashed ray, then 
the heavy ray; but not both. 

The only time the method selects a ray like that refracting at 
D is in the rare instance when point B happens to lie very near the critical 
angle for the lower layer. This feature is generally good because when 
refraction occurs in the vicinity of the critical angle a significant de- 
focusing occurs with attendant loss of signal strength. But if strong 
D-layer absorption occurs on the heavy ray, the defocusing loss at D could 
be less than absorption loss; in this event the method misses the best 
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propagating path. Again, a more elaborate procedure could be devised to 
cure this weakness. 

Other cases of interest which are not handled at present include 
spread F and sporatic E multipath, and cases where unusual spatial depen¬ 
dence of the layers produce freak (but occasionally real) propagation paths. 

In summary, the selected technique is believed to produce the 
strongest propagating mode under most conditions, but has weaknesses in 
that occasionally a potentially significant alternate mode is missed. 

The most important improvement would be to include sporatic E and spread 
F effects. These can be treated to an adequate approximation for many 
cases as a perturbation on the successful ray without incurring a signi¬ 
ficant penalty in computation time or logic. Treatment of other alternate 
rays will incur a noticeable penalty in computation time, and to some 
extent in complexity of logic. 
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3. SINGLE HOP CALCULATION 

After the great circle path for an n-hop mode has been divided 
into n equal segments, each hop is calculated independently, as indicated 
in steps (1) and (2) of Section 2. Here we discuss step (2), calculation 
of a single reflection from a single parabolic layer, and correction for 
refraction in passage through a single parabolic layer. 

3.1 Reflection 

Our technique differs from 
that of Reference 5 substantially 
only in rearrangement of their iter¬ 
ation scheme. For more details the 
interested reader is referred to 
Reference 5. 

The scheme is a prediction- 
correction method based on the quan¬ 
tity 

x = occ cos0 , 

where a is the ratio of broadcast frequency to critical frequency at 
reflecting layer maximum, The factor c is a correction for 

ionospheric curvature, and 0 is defined in the sketch. 

Let h Q be altitude of the layer bottom and a be layer 
thickness, then the apparent reflection height ("equivalent" or "equivalent 
triangular" height), corresponding to an assumed value of x is 
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h' - h 0 + 


ax 0 /l+x\ 

t ln (nr) 


( 2 ) 


Now the tangent to the angle 0 can be calculated from 


tan© = T-r- 


sin<{> 


+ (l-cos4>) 


(3) 


where R is the earth's radius, 
e 


Also, the true maximum ray height can be calculated. 


h = h + a 
o 


(■ - - X 2 ) . 


(4) 


Next the curvature correction is calculated, 




. 2(h -h) „ 2„ 

c = ^ 1 - IFTiTT tan 0 

v e o 


(5) 


Finally, using the trigonometric identity, 


cos0 = (1 + tan^0) 


( 6 ) 


the predicted values of x can be calculated (from (1)) 


To start the iteration we calculate three values of x given 
initial guesses 


x(x x ) , x(x 2 ) , x(x 3 ) 

Assuming the calculated value of x to be a quadratic function of the input 
value, we can solve for two values of input x which should reproduce them¬ 
selves, 
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x = x ± /F /B 
n o v o 


where 

X o = I t X l ‘ x 2 " R(x 3 ' x 2^ / [ (x i " x 2^ " R(x 3 ' x 2^ 
R = [x(x 1 ) - x(x 2 )]/[x(x 3 ) - x(x 2 )] 

B = [xfXj) - x(x 2 )]/[x^ - x 2 - 2 x q (x 1 - x 2 )] 


Fq = - B ( X 1 - X o^ 


If F > x then no solution exists, the ray penetrates the ionosphere, 
o o 


It is necessary to place a logical guard against prediction of 
negative x; if this occurs set 

x = x /2 . ( 8 

n n-1 

It is also necessary to guard against prediction of x^> 1, if 
this occurs set 


*„ ■ (1 - Vi> /2 


We use x : = 0.1, x 2 = 0.5, x 3 = 0.9 and require for convergence 


|x( x n ) - x n l < 10 


So far, no case tried has required as many as ten iterations to 


converge. 
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and can now be recalculated, including refraction effects. 

This procedure is repeated for each layer below the reflecting 

layer. 


19 





4. 


RANGE ADJUSTMENT FOR DIFFERING REFLECTION HEIGHTS 


The method chosen is sufficiently accurate for our purposes but 
is more cumbersome than one would prefer. 

4.1 Method 


For orientation we first work out the case of a flat earth. 



Given the number of hops and an arbitrary set of reflection 
altitudes, h^, we seek the launch angle, 3, and the set of pass ranges, 
r^, which allow 3 to be constant at each reflection point. 

We have 

tan3 = h^/r^ (14) 

Then 


r. = h./tan3 
i i 

and the total range, R, is 

R ■ HH8 L \ 


(15) 


(16) 
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So the launch angle can be found from 


tang = R 1 h. 


and individual pass ranges from (15). 

With the solution to this / 

easy problem in mind, we now proceed 1 Jr 

to search for an approximate solution / 

to the curved earth problem. In this R g / R g 

case all of the surface reflection / 

angles will again have a common value / 

6, as illustrated in the sketch on p 20_ ^ / 

Failing to find a neat closed form 

solution to the problem we proceed to / 

find an approximation which resembles V 

(14) and is accurate for h^/R «1. Such an approximation is 



h [1 - r / (2h.R 1] 

ta " 6 ■ r a » ■.,/£ 


Equation 18 is good to about 10 percent for Ik of 3000 km and 
r^ small compared to R g , and becomes better rapidly as Ik approaches 
realistic reflection altitudes of a few hundred kilometers. 

Proceeding as in the plane case, we solve (18) for r. 


r i R e k^2x^ + (1+x.) 2 tan 2 0 - (1+x 


i )tanB 


where x. = h./R . 

i i e 
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At this point it is clear that a problem exists, for if one sums 
Equation 19 to force agreement with the total range, analogous to the method 
by which (16) was obtained, we shall not be able to readily solve for tang. 

Accordingly, we make further approximations by studying large 
8 and small 8 limits. 

If 8 is large enough (19) can be expanded to obtain 

h. 

r i = (l+h./RJtanir + 0(;h i) * (20) 

Then, summing as before 

1 h ' 

tan8 R L (i+h./R ) ’ (21) 

i i e 

which looks much like the plane case solution. Equation 17. 

This large 8 limit occurs in practice but unfortunately will 
not be exclusive, or even common, in most situations. The small 8 limit 


r. = 2R h. - (R +h.) tan8 + 
l \/ e i e i 


(R +h.) 2 tan 2 8 
e i _ 

2\/2R h. 

” e i 


Summing (22) to find total range and solving for tan8 yields 


tan8 = C/B -J C 2 /B 2 + 2(R-A)/B 
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where 


i S 

B E £ (R +h.) 2 /\l2R h. / (24) 

. e 1 y e i I 

i 1 

C = £ (W ) 


Behavior of the two limiting solutions can be conveniently 
studied by choosing a case where all h^ are equal, then the large g 
limit. Equation 21, becomes 


tang = 


Nh 

R(l+h/Rj 


(25) 


where N is twice the number of hops. 


The small g limit. Equation 23 becomes 


tang = 


V2Fh 

( R e+ h) 


1 - 


2R 


Define 


E -\j 2N 2 hR e /R 2 


N\2Rh 


(26) 


and 


(27) 


F E N(R e +h) tang/R 
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then Equation 25 becomes 


F l = q 2 /2 (28) 

and Equation 26 becomes 

F s = q (1 ->/2/q-l) (29) 

The behavior of F vs. q is illustrated in Figure 2, the exact 
solution illustrated for this special case is obtained readily from Equation 21. 

A simple technique with accuracy adequate for our purposes, 

suggested by Figure 2, is to calculate both limits, that is, Equations 21 

and 23, then choose the smaller of the two. 

For realistic cases, q is never as small as 0.1 (which implies 
h < 1 km) so the error near q = 0 is not significant. The maximum 
fractional error of about 30 percent occurs at q = 2, typically at the 
very high launch angle of 30°. 

Much of the error is removed by a final adjustment which forces 
the r^ to precisely sum to the correct range. 

4.2 Example 


An example which stresses the algorithm somewhat is summarized 
in Table 1. The transmitter is located at Washington, D.C. and the re¬ 
ceiver at Tel Aviv, at 84.92° great circle distance (9432 km). Ionospheric 
conditions are appropriate for 0600 UT, which places the western two-thirds 
of the link in darkness (local time in Washington is about 1:00 a.m. and 
in Tel Aviv is about 8:30 a.m.) 
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Table 1. Height adjustment for a four-hop mode. 

Column 1 lists the pass, or break point index (1 is 
Washington), columns 2 and 3 list north latitude and 
east longitude of equally spaced break points. The 
fourth column shows altitudes, followed by break point 
coordinates adjusted by the algorithm to account for 
differing h r . The seventh column lists ground angles 
which correspond to the adjusted coordinates, followed 
by corresponding pass range and the error in break point 
position due to our approximation. All angles are in 
degrees and all distances in km. 


1 Break 
Point 

Equal Spacing 

r 

Adjusted for 

hr 

Pass 

Range 

Position' 

Error 

Lat. 

Long. 

h r 


Lat. 

Long. 

Angle 

1 

33.90 

-77.03 

0 


38.90 

-77.03 

6.97 

1260.81 

0 

2 

44.84 

-65.19 

290.30 


45.20 

-64.28 


1260.81 

1.88 

3 

49.30 

-51.01 

0 


49.77 

-48.86 

6.96 

1310.66 

3.76 

4 

51.74 

-34.72 

307.53 


51.98 

-30.43 


1310.66 

6.43 

5 

51.75 

-17.54 

0 


51.15 

-11.44 

6.92 

1598.57 

9.10 

6 

49.35 

-1.23 

417.43 

j 46.34 

9.22 


1598.57 

6.59 

7 

44.91 

13.00 

0 

! 38.51 

25.65 

7.05 

548.21 

4.08 

8 

38.99 

24.86 

92.67 

j 35.37 

30.39 


548.21 

2.04 

9 

32.06 

_ 

34.76 
- ■ ■ ■ 

0 

; 32.06 

34.76 



0 
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The mode illustrated has 4-hops, three F2-layer reflections 
and one E-layer reflection. As the ray traverses the great circle path 
successive reflections occur at increasing effective altitudes until the 
last reflection, which is at sharply decreased altitude. The reasons for 
this behavior illustrate the complicated interplay of ionospheric behavior 
and radio propagation. The second reflection point is above the first 
predominantly because the F2 and FI layers are weaker at the second hop 

(break point 4 in Table 1) where the local time is about 4:15 a.m. The 

third hop (break point 6) has higher effective altitude despite the dawn 
buildup of the F2 layer, because the FI layer is nearly strong enough to 
reflect the ray and thus refracts strongly. Finally at the fourth hop 
(break point 8) not only is the FI layer capable of reflecting the ray 
but also the E-layer can reflect and our algorithm chooses the lower mode. 

An exact solution would have yielded equal angles approximately 
7.008°. Thus the seventh column of the table can be used as a measure of 
inaccuracy of our algorithm. Another measure, one which is more directly 
pertinent to nuclear effects, is the error in placement of the break points. 
This error measures the distance by which a disturbance, beta patch or 

other, might be missed and is shown in the final column. The error is 

seen to build up to about 9 km then fall back toward zero, despite the 
fact that break points were shifted as much as 1397 km (No. 7). Since 
any natural or nuclear disturbance of interest will be large compared to 
9 km one concludes that the approximation is adequate. 
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5. 


AMBIENT IONOSPHERIC TILT ANGLE CALCULATION 


The main motivation for inclusion of ambient ionospheric tilt 
effects is to establish a mechanism by which tilt effects in a highly dis¬ 
turbed ionosphere can be calculated. For this reason we neglect lower 
refracting layer effects and treat only reflecting layer tilt. 


In the sketch the solid 
lines illustrate real ray paths, 
one of which connects points (1) 
and (2). This ray reaches a max¬ 
imum height, h m> and is to be 
represented by an equivalent 
triangular path which reflects 
at height h r> illustrated by 
dotted lines. To find the effec¬ 
tive tilt angle, t, of the 
reflection point, we calculate 



the true maximum height, h', of a similar ray, with points (1) and (2) 
displaced a distance 6r. 


m 
Then 


T = 


(h -m 

m m 
6r 


(30) 


Equation 30 is good to first order which is adequate for the first 
treatment of this effect; we do not hope for accuracy when t is large. 

For the moment, and perhaps forever, it will have to be sufficient to know 
that the mode is highly disturbed when T is large, and the approximate 
degree of disturbance. 
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5.1 Tilt Parallel to the Great Circle Path 

After calculation of the first approximation to reflection, de¬ 
scribed in Section 3, the true maximum ray penetration height, h^, is 
known from Equation 4 and the corresponding value of x from Equation 1. 

It only remains to choose one or more other points along the great circle 
path to calculate h^, using the same value of x. Since the ionosphere 
is described by parameters specified at points along the great circle, the 
technique adopted is to locate the two great circle ionospheric data points 
surrounding the reflection point, evaluate h m at these two points, then 
calculate "parallel tilt", T„ from (30). If the ray is capable of pene¬ 
trating the ionosphere at one of the great circle data points, the reflec¬ 
tion point is used in its place. The convention is adopted that tilt is 
positive if height increases toward the receiver. 

5.2 Tilt Perpendicular to the Great Circle Path 


It would seem that our method of ionospheric description, which 
restricts knowledge of ionospheric parameters to those along the great 
circle path, might prevent finding tilt angles perpendicular to the great 
great circle. This is not the case because at each point on the great circle 
path the ionospheric data table covers all 24 hours of local time. This is 
equivalent to data at all values of longitude provided effects of geomag¬ 
netic coordinate changes on the longitude dependence in the neighborhood 
of the point of interest are neglected. 


The situation is depicted ^2 
in the sketch. We know the latitude, 
longitude coordinate of the reflec- X 
tion point (X,Jl), and the two great 
circle ionospheric data points 
surrounding it (A^,^) and ^1 

(A 2 ,)l 2 ). The tilt angle 
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p 


perpendicular to the great circle path can be obtained by differencing 

h at the points (A,,£ ) and (A„,£ ). 
m r 1 e 2 p 

Since arc lengths between great circle data points are small, plane 
geometry will suffice. In this approximation the interesting longitudes are 
given by 


= £ + (A-A 1 )(A 2 -A 1 )/(£ 2 -£ 1 ) 

(31) 

= £ + (A-A 2 )(A 2 -A 1 )/(£ 2 -£ 1 ) 

(32) 


If the local time at (X,H) is t, then ionospheric parameters 
for (Aj,£ e ) can be obtained by evaluating those at (A^,£^) at 

t = t + a(£ g -£) (33) 

where a is the factor appropriate to convert longitude to time (in our 
case ionospheric time is measured in radians from local midnight, so a = 1). 
Similarly 


t = t + a(£ -£) 
P P 


(34) 


is the local time appropriate for point 


(vy 


Then the perpendicular tilt, t^, is obtained from 


T i = t h m O) “ h m CpDl/ d ( 35 ) 

where d is the distance between points (^j»*- e ) and ^2’^p"*' Note 
convention adopted here is that perpendicular tilt is positive if height 
increases to the right when looking from transmitter toward receiver. 
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There are two special cases. If the path segment 1 ■+ 2 is very 

nearly North-South < 10 rad) we evaluate h^ at (X,£ + 20 

minutes) and (A,£ - 20 minutes). If the path segment 1-► 2 lies too 

2 

nearly East-West (IX^-X^I < 10 rad) this scheme fails and zero perpen¬ 
dicular tilt is assumed. 
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6 . 


EFFECT OF IONOSPHERIC TILT ALONG THE GREAT CIRCLE PATH 


Consider the case of a propagating mode with arbitrary but fixed 
number of hops where the ionosphere varies along the path such that each 
hop has a different virtual reflection height, h^, and has a small tilt 
angle, x^; but the ray remains in the great circle. As before, the path 
composed of straight line segments which are called "passes" connected at 
"break points". Thus an n-hop mode consists of 2n passes and 2n + 1 break 
points as depicted in the sketch below, for a 4-hop mode. 



Transmitter 


Receiver 


The tilt angles are given by 



(36) 


The entire mode geometry has been calculated as previously de¬ 
scribed including effects of varying reflection height, but under the 
assumption that the ionosphere is circularly symmetric about the center 
of the earth at each reflection point, that is, effects of tilt have been 
neglected. A formalism is needed to modify break point coordinates to ac¬ 
count for tilt, one which is accurate for reasonably small deflections and 
"sensible" for large. By "sensible" we mean the break point range should 
be mom tone increasing with break point, reflection altitude shouldn't 
suffer a major fractional change, etc. 







6.1 


Single Hop - Plane Geometry 


Initially consider the problem of a single hop in plane geometry, 
then one can work toward a multihop problem in circular geometry. In the 
sketch below, dashed lines show the situation which would obtain if the 
ionosphere were parallel to the surface of the earth while solid lines 
show the situation which obtains when the ionosphere has a (positive) tilt 
angle T. 



From the sketch, the launch angle without tilt is 

ctng = tan a = ^v- (37) 

o 2h ' 

and, if launched at angle 6 , the ray will be land beyond the receiver at 
angle 6-2x. 
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ctn (8-20 - i * -jf 


(38) 


As a first step in finding the adjustment in £, h', and r 
required to force the ray to reach the receiver, search for an expression 
for the angle 


a = (Yj + Y 2 )/2 
as a function of tilt angle t. 

Bisection of y by the dashed line shows 
Yj = a - t and y 2 = a + x 
One useful relation can be found by combining 
h'tan (a-T) = D/2 - 6 

and 

h'tan (a-t) = D/2 + 6 

as 


(39) 


(40) 


(41) 


(42) 


h'[tan (a+x) + tan (a-x)] = D . (43) 

In order to eliminate h' another relation is required. Use 

tanx = (h-h')/6 ( 44 ) 

with (41) to eliminate <S, 
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(h-h') ctnx = j - h' tan (a-x) (45) 

Elimination of h between (45) and (43) followed by an unusually tedious 
amount of trigonometric manipulation leads to the startling result 

h tan a = D/2 ! ( 46 ) 

In short, the angle by which the successful ray is reflected by the tilted 
ionosphere is equal to the untilted value, independent of X. At risk of 
belaboring the obvious we point out that this result is not a small angle 
approximation but is exact. 


Having established that the tilted and untilted cases have a 
common apex angle it is a short exercise in plane geometry to show that 
the launch angle has been raised 


SB t = x 


(47) 


while the receiving angle has been decreased by x. 

From (43), now that a is known to be constant, find the change 
in reflection altitude. 


2h tana 


h'- h = 


tan(a +x) + tan(a -x) 
o o 


- h 


(48) 


After more exercise in trigonometric substitution. Equation 48 simplifies to 



(49) 
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Using (44) one finds the ground range of the reflection point to 
have been shifted by an amount 


h-h' 

tanr 


sint cost 
sin 2 3 


(50) 


where the minus sign is adopted to show that positive X (reflection alti¬ 
tude increasing with ground range) produces a shift toward the transmitter. 


This single hop case will be useful primarily as a check on approx¬ 
imate methods to be developed for more complex cases. Attempts at preci¬ 
sion for more complex cases rapidly become bogged down in more trigonometry 
than justified on the basis of our assumption of constant x in the 
vicinity of a given break point. 

6.2 Heights Constant - One Place Tilted - Plane Geometry 



Now consider the cumulative effects of multiple hops on break 
point range. In the sketch a ray has been launched at angle G appro¬ 
priate for a multiple hop path to the receiver at (R) , the n t *‘ break 
point. 
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The tilt T causes reflection angles to change by 2 t after 
point m. Thus the distance, 6r, by which the receiver is missed is 

6r = (n-m) h(ctn(B-2T) - ctnB) 

** (n-m) h(l+ctn 2 6)2T . ( 51 ) 

The total range from (t) to (R) is 

R = (n-1) h ctnB . (52) 


then 


HR ? 

gg- = - (n-l)h(l+ctn B) . (53) 

So the change in launch angle required to bring the ray down on the receiver 
is found from requiring 

§ 66 - - 6r 


and is 


dB 


2 (n-m) 
(n-1) 


(54) 


We can check our approximation by noting that for a one-hop case, 
n = 3, m = 2, (54) gives 


dB = t , the exact answer, while for an extreme case, m 

small, n 


dB -* 2t, 
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Reference to the sketch shows the angle 2 t to be the one which yields 
the unperturbed angle 3 thereafter. For this extreme case the approxi¬ 
mation has introduced an error equivalent to the miss distance on one hop, 
an error commensurate with current understanding of the ionosphere, especi¬ 
ally so because some of the error can be eliminated by retroactive adjust¬ 
ment of break points to force the ray to land on the receiver. 

The increment in range at intermediate break points is given by 
/h ctn [3 + 2 T] , i^m 

<5 r i = \ 

(h ctn [3-2 T 1 > i > m (55) 


6.3 Heights Vary - One Place Tilted - Plane Geometry 


To find the miss distance we must sum heights in Equation 51 
n 

6r = h. (l + ctn 2 3) 2 t 
n m 1 

n 

= 2 (l+ctn 2 B)T £ Ih -h I . (56) 

m+1 1 1 

Following an argument similar to the previous case we find 



(57) 
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6.4 


Heights Vary - All Places Tilted - Plane Geometry 


From (57) the change in launch angle due to tilt at location 


m is 


60 = T 

m m 


n 

2 £ hi 
i=m+l 
n 

£ \ 

i=2 


(58) 


The total change in launch angle is given by summing (58); 

n 


60 = 


n-1 

2 £ 2>i 

m=2 i=m+l 


n 

m= 2 1 


(59) 


The simplest method of finding the new ranges and heights of 
break points is to accumulate them, viz. 
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(6 


and 



fr'-r )tan x 
mm 


m 


where a second order effect in the equation for r has been neglected. 

m 


When one reaches the last break point the accumulated error is 
measured by the miss distance. All ranges are then scaled to force agree 
ment in total range. 

6.5 Extension to Circular Geometry 

If 3 is interpreted as the launch angle relative to the chord 
of the arc as depicted below nothing changes to first order in the pass 
range over the radius of the earth. 



US 





It is necessary, however, to rederive (59) accounting for the fact 

that, while incident and reflection angles are equal relative to the earth's 

surface, they are not equal relative to the chord of the arc. The effect of 

2 

this extension is that the factor (1+ctn 3) cannot be regarded as constant, 
but must be retained. Equation 59 becomes 


63 


n-1 
2 £ 
m=2 


T m X Ml + ctn 


i=m+l 




£ h.(l + ctn 2 3 i ) 
i=2 


(62) 


In many cases a second order term due to the change in incident 
angle which results from change in pass range is significant. Proper treat¬ 
ment of this term is a bit cumbersome and more importantly, tends to de¬ 
stabilize the adopted mathematical treatment. A satisfactory expedient is 
to retroactively adjust all pass ranges, to force agreement with total range, 
by an amount proportional to the square of the calculated pass range. 


6.6 Example 


Table 2 summarizes results of the approximate treatment developed 
above when applied to the same four-hop Washington to Tel Aviv mode illustrated 
previously. This particular mode was chosen primarily because it represents 
as difficult a case for the parallel tilt approximation as has been encountered 
to date. The table shows errors as large as about 20 km in placement of 
break points. Since this error is small compared to the size of any possible 
D-layer absorbing region the method is considered satisfatory. The calcu¬ 
lated launch angle at Washington is 5.11° compared to the exact solution of 
5.01° and the calculated received angle at Tel Aviv is 12.15° compared to 
the exact solution of 11.89°. 


41 









Table 2. Parallel tilt adjustment for a 4-hop mode. 


Column 1 lists the break point index (1 is 
Washington, D.C., 9 is Tel Aviv). The next 
three columns list latitude, longitude, and 
great circle range measured from Washington 
after adjustment for different reflection 
height. The following four columns show the 
tilt angle parallel to the great circle path, 
reflection height and position adjustments 
for tilt and the increase in range due to tilt. 
The next column shows the exact trigometric 
position (accepting h~) and the last column 
the amount by which tne approximate method 
overestimates the position. All distances 
are in km and all angles in degrees. 

















7. TILT ACROSS GREAT CIRCLE 

7.1 Plane Case 

Looking from transmitter toward receiver, one sees a tilted ionos¬ 
phere and the ray path. Positive tilt means h increases to the right as 
depicted below. Launch a ray along the great circle path. Without tilt 

one would see only one vertical line. With tilt n at break point m 

m r 

the ray is deflected to the right. 



The ray comes off the ionosphere at angle 2n , causing it to land a distance 

m 

21^ to the right of the great circle path, then if there are no more 

tilted break points, it moves additional distance 2p h. for each pass. 

m j r 


The total displacement at the end of the path, break point n, 
due to the tilt at break point m is 


d 

m 


2n 


m 


n 


£ 

i=m+l 


I h. -h 



(63) 


If one assumes a tilt at each break point, one has 


n-1 

n-1 

n 


£ ■ 
m=2 

2 £ n m 
m=2 

V |h.-h. ,1 

i=m+l' 1 

(64) 
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The launch azimuth (positive to right) of the ray which reaches 
the receiver must counteract the deflection; 


a L - - d/R 


where R is the great circle range of the total path. 

The azimuth of the ray at the receiver is 


a = - d . / r 
R n-1 n 


where r^ is the ground range of the final pass. 


Since the launch azimuth is a., the resulting displacement 
perpendicular to the great circle of the J break point is 


J-l 

J 


s T = 2 Y n 

Y h.-h. + ad 

.*-> l l-l L J 

(67) 

m=2 

i=m+l 1 1 


where d. is ground range to 

point J. 



The new reflection altitude is 


h J = h J + S J tan n j 


n J-l 

Numerically it is convenient to collect Y |h.-h. ,1 and Y P , then 

i=m+l 1 1 m=2 

s, * s j.i* ">j-i h i h i- h i-ii * 2 i h J- h J-iiM n » * “ l rj 




i -h T , V n ♦ a. r 
J J-l m L J 


where r^ is the incremental range of the J pass. 





7.2 Extension to Spherical Geometry 

If Ik is interpreted as the height of reflection above the chord 
of the hop rather than above the earth's surface. Equation 69 is still valid. 

7.3 Determination of Break Point Horizontal Coordinates 



If A is colatitude and £ is longitude, to a good approxima¬ 
tion one can write 


Sj = R e ^(<5A) 2 + sin 2 Aj (6£)‘ 


(70) 


Since the slope of the line segment Sj is the negative gradient 


of the line segment r T , 

J 


dA 

d£ 


■ - (§) - - 
'T 


- 1 'V'j-i) 




(71) 


combining (70) and (71), 

6A = 


± S, 




(72) 
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and. 




where the positive sign is chosen if £ > £ , 

J “ I 


± S, 


<5£. = 


J * si " 2x j 

where the positive sign is chosen if X < X. 


J-l ’ 


(73) 


7.4 Example 

Table 3 illustrates results obtained for the Washington-Tel Aviv 
case. The distance by which the ray missed the receiver (1.5 parts in 10 ) 
measures the combination of second order error in the technique and round 
off error in a 32 bit computer word length. Of more interest is the much 
lower magnitude of perpendicular displacements relative to parallel dis¬ 
placements. Most of this reduction in effect is due to the circumstance 
that the path lies mostly East-West so that perpendicular gradients are 
smaller than parallel gradients and the further accident that perpendicular 
tilts tend to cancel (the first two angles sum to +0.26 and the next is 
-0.26). However, parallel gradients systematically tend to be about an 
order of magnitude more effective in displacement than perpendicular gra¬ 
dients. The parallel effect tends to be enhanced by small angle of incidence 
at the earth, while the perpendicular effect always acts as though the ray 
were near a right incident angle. 
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Table 3. Perpendicular tilt adjustment for a 4-hop Mode. 

Column 1 lists the mode break point index, column 2 
reflection altitude, and column 3 the tilt angle 
perpendicular to the great circle path. Column 4 
lists displacement of break point, positive to right 
of great circle path as viewed from the transmitter, 
and column 5 the adjusted reflection height due to 
perpendicular tilt. All distances are in kilometers 
and all angles in degrees. 



h r 

- - "" - ' i 

n 

After Adjustment 

s 

h 

i 

0 

0 

0 

0 

2 

290.03 

0.11 

-1.29 

290.03 

3 

0 

0 

-1.40 

0 

4 

308.50 

0.15 

-1.49 

308.49 

5 

0 

0 

0.02 

0 

6 

408.41 

-0.26 

2.34 

408.40 

7 

0 

0 

0.91 

0 

8 

92.54 

-0.03 

0.43 

92.54 

9 

0 

0 

-0.14 

0 
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